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ABSTRACT 



We use the Dirac equation with a "(asymptotically free) Coulomb + (Lorentz 
scalar) linear " potential to estimate the light quark wavefunction for qQ 
mesons in the limit rriQ — ► oo. We use these wavefunctions to calculate the 
Isgur-Wise function £(v.v') for orbital and radial ground states in the phe- 
nomeno logically interesting range 1 < v.v' < 4. We find a simple expression 
for the zero-recoil slope, £'(1) = — 1/2 — e 2 < r 2 > /3, where e is the energy 
eigenvalue of the light quark, which can be identified with the A parameter 
of the Heavy Quark Effective Theory. This result implies an upper bound of 
— 1/2 for the slope £'(1). Also, because for a very light quark q(q = u, d) the 



size J < r q 2 > of the meson is determined mainly by the "confining" term in 
the potential (7 D crr), the shape of £ u ,d(v-v') is seen to be mostly sensitive to 
the dimensionless ratio A 2 J a. We present results for the ranges of param- 
eters 150MeV < K u4 < 600MeV (A s « A u4 + lOOMeV), O.UGeV 2 <a< 
0.2hGeV 2 and light quark masses m u , fa 0, m s = YlhMeV and compare to 
existing experimental data and other theoretical estimates. Fits to the data 



give: Al >d /a = 4.8±1.7, -£^(1) = 2.4±0.7 and 1^1^^ = 0.050±0.008 

[ARGUS '93]; k\J a = 3.4 ± 1.8, -^(1) = 1.8 ± 0.7 and \V*\^E: = 

0.043 ± 0.008 [CLEO '93]; A 2 u4 /a = 1.9 ± 0.7, -£' u>d {l) = 1.3 ± 0.3 and 
\V cb \F(l) = 0.037 ± 0.003 [CLEO '94] (Existing theoretical estimates for 
F(l) fall in the range 0.86 < F(l) < 1.01). Our model seems to favour 



the CLEO '94 data set in two respects: the fits are better and the resulting 
ranges for the model parameters (A Uj(i , a) are more in line with independent 
theoretical estimates. 
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1 Introduction and Summary 



During recent years, a lot of effort [Ij has gone into the description of systems and pro- 
cesses involving at least one heavy quark (juq ^> Aqcd) via a systematic expansion in 
the small parameters (Aq C D/ m Q) and a s (mg), taking advantage of "Heavy Quark Sym- 
metries". The usefulness of this approach is in that in the heavy quark limit niq — > oo, 
all the physics can be expressed in terms of a small number of form factors, which depend 
on the light quark and gluon dynamics only. These "universal" functions can be used 
as a means for comparison among different theoretical models (such as non-relativistic 
and relativistic potential models, lattice QCD calculations, etc). Comparisons with ex- 
periment or with methods such as QCD sum rules require in general introducing model 
dependent O ( ) corrections ( except at v.v' = 1, where corrections start at order 

(Aqcd 2 l m Q) 0)> as wen as calculable perturbative QCD corrections. 

In this paper we limit ourselves to the leading order in the heavy quark expansion. 
We calculate the Isgur-Wise function ||, £(v.v'), for radial and orbital ground state 
mesons using a relativistic model for the light quark. This function in our case (radial 
and orbital ground state), for a given flavour of light quark, fully describes the transition 
M q — > M' q in which a local operator transforms the heavy quark Q in the initial meson 
M q (of 4-velocity v ; J=0 or 1) into the heavy quark Q' in the meson M' (of 4 velocity 
v' ; J=0 or 1). 

We assume that the light quark wavefunction obeys a Dirac equation with a spheri- 
cally symmetric potential in the reference frame in which the heavy quark is stationary at 
the origin. We also assume that the potential has the form V(|x|) = K(|x|) +c +7°cr|^|, 
where V^(|x|) is an asymptotically free Coulomb term, c is a constant and the last term 
is a Lorentz scalar confining term. Thus, the spatial wavefunction for the light quark q, 
of mass m q , is assumed to obey the time-independent Dirac equation 

a-HVa?) + V c (\x\) + c + 7>|£| + m q )] = e*(£) . (1) 

Because we are investigating the meson system in the heavy quark limit (AqcD/ m Q 
0), the energy eigenvalue e of the Dirac equation can be identified with the "inertia" 
parameter A q often introduced in the Heavy Quark Effective Theory [HQET] 

In our investigation, we allow the A parameters to vary over the range 0.15 < A u ^ < 
O.QGeV (A s ~ A Uj d + lOOMeV), obtained from recent lattice gauge theory studies f| 
and other theoretical estimates 0. As will be seen below, the parameter A plays 
an important role in the calculation of the Isgur-Wise function. The additive constant 
c in the potential allows us to pick an arbitrary value for A. The asymptotically free 
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Coulomb term V c (r) is parametrized in terms of the QCD scale Ajjs an d a saturation 
value for the strong coupling af = a s (r = oo) (We use Ajj$ = 240MeV and af ~ 1). 
Phenomenological and theoretical considerations motivate our use of 3 different values 
for the string tension : a = 0.25GeV 2 , Q.18GeV 2 and O.lAGeV 2 . We set m u , rrid = and 
m s = 0.175GeV . Further details about the parametrization of the potential are given 
later in the paper. 

We find that the shape of the Isgur-Wise function is mostly sensitive to the parameters 
A and a. In fact, to a very good approximation £, u ,d(v -v') depends on the dimensionless 
ratio — only Within our approximation, we also find a strict upper bound on the 
slope at zero recoil: £'(1) < —1/2. A recent paper |7| that uses the MIT Bag Model 
formalism obtains the same bound, which is stronger than the well known Bjorken bound 
£'(1) < —1/4 |5|. For large values of v.v' (2 < v.v' < 4) the shape of £(v.v') is relatively 
insensitive to the input parameters A, o and Ajjg. 

We present explicit results for £(v.v') in the phenomenologically interesting region 
1 < v.v' < 4, using the above mentioned range of parameters. We compare our results 



with recent experimental data [g, [T(J, [TT| of semileptonic B decays as well as with other 

A 2 

theoretical estimates. From fits to the ARGUS '93 data [0] we obtain the ranges — — = 



4.8 ± 1.7 (corresponding to -£^(1) = 2.4 ± 0.7) and Jj^- S = 0.050 ± 0.008. 
Although the fits are of good quality, for reasonable values of the string tension a 
they favour values for the inertia parameter A u ^ that are significantly above the range 
Q.lhGeV < A Ut d < 0.6GeV advocated by most theoretical estimates [|], |, |j. The 
corresponding range for the slope — £u,d(l) is also above most independent theoretical 
estimates (see discussion in section 4). The best fit to the CLEO '93 data |10| was of 

A 2 

significantly poorer quality. Here we found the ranges = 3.4 ±1.8 (corresponding 

A. 2 



to -&,«*(!) = 1-8 ± 0.7) and \V cb \^ j^- g = 0.043 ± 0.008. The ranges for has in 
this case some overlap with previous theoretical estimates, but is still somewhat on the 
high side of these estimates. We noticed that if the data point from the CLEO '93 set 
[|T0| corresponding to the largest value of v.v' is ignored, the quality of the fit greatly 
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improves. The parameter ranges obtained in this case are = 2.1 ±1.6 (corresponding 
to -£^(1) = 1.3 ± 0.6) and iKblyofb; = 0.038 ± 0.006. These ranges for ^ and 
— £4<z(l) overlap with many previous theoretical estimates, but are somewhat too wide 
to provide useful new information that could distinguish between these. 

The recently released CLEO '94 data analysis ]TTJ] has smaller error bars than the 
data mentioned above 0, |TU|| . Also, our model produces better quality fits to this 
data set than to the previous ones. These two factors contribute to help narrow the 

A 2 

parameter ranges. The following ranges are favoured: = 1.9 ± 0.7 (corresponding 
to -Cd( l ) = 1-3 ± 0.3) and |U c6 |F(l) = 0.037 ± 0.003 (Different theoretical estimates 
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for the constant F(l) can be found in eq (40), section 4). Note that the central values 
of these ranges are quite close to the ones found from the fits to the CLEO '93 data 
set after removing the CLEO '93 data point corresponding to the highest v.v'. The 
uncertainties, however, are here reduced by a factor of about 2, which should help us to 
narrow down the parameter space of the underlying physics. 



2 General Formalism 

We assume that the wavefunction \1/ of a qQ meson of mass M q q in the limit uiq — > oo 
(i.e ignoring 0(^2^) effects) can be expressed in terms of a direct product of the (free 
spinor) wavefunction of the heavy anti-quark, x, and the wavefunction of the light quark, 
■0, in the "relative" coordinates. It is therefore convenient to define the 4-vectors 

X» EE {t Q , f Q ) , X» EE (t q - t Q , f q - f Q ) (3) 

Of course, the kinematics of the two wavefunctions x an d i> are not completely indepen- 
dent. They are connected by the implicit constraint that in the rest frame of the heavy 
quark (which is also the rest frame of the qQ meson) the spatial part of the light quark 
wavefunction, ip(x), obeys the time-independent Dirac equation (1). In order to be able 
to give the meson wavefunction \1> a physical meaning in a given reference frame, we 
have to set the "relative time" of the two constituents to zero 



X° = , t q = t Q = t 

In the rest frame of the meson (P = 0), its wavefunction can then be written as: 



2M, 



qQ 



\i w 



A C 



-iet a 



(4) 



where A , 77 are spin indices (A , r] —\ or \ with respect to some axis); xjj is a f ree 
Dirac spinor at rest corresponding to the heavy antiquark and ip^{x) obeys equation (1) 
with energy eigenvalue e. The factor in front of the direct product is for normalization 
purposes (See eq (10) below). Using our convention tQ — t q — t we see that the rest 
mass M q Q is given by 

M gQ = m Q + e (5) 

thus validating our identification of e with the parameter A q commonly introduced in 
HQET. 

We now turn to the description of the meson qQ as seen from a reference frame L, with 
respect to which the meson (rest frame L') is moving with a 4- velocity v = (7, 7/9). In the 
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unprimed frame, the 4-momentum of the meson is given by P = vM. We assume that 
the wavefunctions of the light quark and heavy antiquark transform in the standard way 
under a Lorentz transformation. Although this would be strictly correct only if the Dirac 
equation for the light quark (1) were in fact covariant under Lorentz transformations, 
we assume (and verify in some special cases later) that this procedure will not introduce 
important errors in our results. We can then express our wavefunctions in the unprimed 
frame in terms of those in the rest frame of the meson, V 



Xv(tQ,r Q ) = S v Xd(t' Q ,r' Q ) 



(6) 



where 



t q ,Q - P-r q ,c 



7„ 



f 9i Q + ( 7 -l)(/?.f 9 ,o)/3-7/3t 



q,Q 



tn=t 



Q 



and 



S v = exp 



tanh~ x \(3\ 



h + i 



1 








1 






7+1 


7+1 


y(/3 x +i/3y) 


7/3* 


7+1 


7+1 



7+1 

7(Ae+»A/) 
7+1 

1 





y(j3x-i!3y) 

7+1 
_7A 
7+1 





(7) 



Therefore, the wavefunction of a meson moving with velocity (3 in our (unprimed) frame 
will be 



2AI. 



2M qQ r -iP.X 



^vXn 



(9) 



where we have used eq (7) for the last step and t q = tQ = t, fq = X , f q — fq = x , P = 
( 7 M, 7/3M). 

It is straightforward to check that the standard normalization for meson states, 

< >= 2P°5\P - P')6 XX >6 W , (10) 



is obtained provided that the usual normalization is used for the light and heavy quark 
wavefunctions: 



Xq Xq 



5 x\' and / d s xip^ (x)ip^ (x 



(11) 



The Isgur-Wise function £(v.v') (for orbital and radial ground states) can be ex- 
tracted in a simple way [JTJ, ^ by calculating the following matrix element between two 
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pseudoscalar states of mass M and 4-velocities v and v'\ 



t{v-v') = {l + vv , )M < PvK(0)h v (0)\P v > , (12) 



where \P V >= ^ > -|*i T >J can be obtained from eq (9) and h v >(0)(h v (0)) is 
the creation (annihilation) operator for a heavy quark of 4- velocity v'(v ) at the origin of 
space-time and v^j^hy = h v has been used. The light quark wavefunctions ip 1 ! (see eq 
(9)) correspond in this case to the radial and orbital ground state solutions of the Dirac 
equation (1). 

Without loss of generality we can assume that (3 and (3' are collinear. We can then 
use the simple result 

Xq ] S v > ] i°S v Xq = J 1 + ^' V (13) 
to obtain the explicit expression 



= ^Jy^j J d'x\ £ ^\x')S^S v ^')e^^'^ , (14) 

where 



x" = x+^-l)((3-x)(3 



(15) 



and v.v' = 77' [1 — (3.(3'] (v.v' in [1, 00]). x 1 and x" represent the spatial coordinates in 
the reference frames where the mesons of respective 4-velocities v' and v are at rest. 

In principle, the above integral expression (14) should depend on v.v' only, i.e. should 
be Lorentz invariant so that its value be independent from the (unprimed) frame that is 
chosen to evaluate the integral in. In practice, because the Dirac eq (1) is not Lorentz 
covariant (due to the potential), the wavefunctions are not either. Therefore, the function 
£(v.v') does depend on the frame we choose. However, it can be easily checked that the 
important properties £(1) = 1 and Im(£,(v.v')) = are satisfied in any Lorentz frame. 
Moreover, as will be detailed at the end of this section, we found that the value of the 
slope at zero recoil £'(1) i s the same in three simple (but quite different) Lorentz frames. 
Thus, at least in the vicinity of the zero recoil point, the effects of Lorentz non-invariance 
of our function £ are expected to be small. 

We choose to work in the Breit [0] frame, where it is easiest to perform the calculation. 
In this frame the incoming and outgoing mesons are moving with equal speeds but in 
opposite directions {(3 = —(3') and therefore v.v' = 2^ 2 — 1. Also, S v , = S~ x so that after 



a change in integration variable £(v.v') can be written in the simple form [Jl~2| , with (3 
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and /3' in the z direction, 



7 

where 



Uv-v') = - 2 J £ \wi^ Z ' > ( 16 ) 

7 T?=Ti 




7 = V— • /3 = V^n- (17) 

2 

Because X^=u ^(^O is spherically symmetric, the angular integration is trivial. 
Also, it is easy to derive an expression for the slope at zero recoil £'(1) in terms of the 
inertia parameter e q = A q and < r 2 q >^: 

= -\-y>,<r\>s=-\-\^<r\>- a . (18) 

This is a general result within our formalism and not linked to any particular form of 
the potential used in the Dirac equation for the light quark. It tells us that £'(1) depends 
only on the light quark energy eigenvalue e (or equivalently for our model, the HQET 
"inertia" parameter A q = lim mg _ >00 (MQg — ttlq)) and on the rms distance between the 



light quark and the (stationary) heavy antiquark, ^J< >^. We should point out that 

this same result was obtained independently in a recent paper that describes qQ mesons 
and qqQ baryons in the context of the MIT Bag Model JF), also obtained in the Breit 
frame. 

An interesting aspect of this result is that it sets an upper bound on the slope (£'(1) = 

V) 

£'(1) < -\ (19) 

This bound is stronger than the well known Bjorken bound £'(1) < — 1/4 and is a result 
of the dynamical assumptions inherent in our treatment of the meson system in the 
heavy quark limit. In particular, we think that it can be traced to the fact that we 
describe a moving meson by boosting the light and heavy quark independently (by the 
same amount). 

Because our formalism is not fully covariant (spherically symmetric Dirac equation 
potential is put in by hand) one may suspect that the above result for £'(1) ( ec L 18) i s 
Lorentz frame dependent. We have checked explicitly that the result is unchanged if 
£'(1) is calculated in the frames where either the incoming or the outgoing meson is at 
rest. 
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3 Parametrization of the Dirac equation potential 



In the Mq — > oo limit, the qQ meson system should be well described by the heavy anti- 
quark stationary at the origin and the light quark (of mass m q ) moving in a spherically 
symmetric static external potential. We are of course ignoring the self interactions of the 
light quark in the hope that these can be absorbed to some extent in a renormalization 
of the parameters of the external potential. 

We take the short distance behaviour of the potential from renormalization-group- 
improved QCD perturbation theory and the long distance behaviour from lattice and 
other non-perturbative studies that ignore screening by light quark pair creation. Unfor- 
tunately, knowledge about the leading behaviour of the potential in these two extreme 
distance regimes defines the potential only up to an additive constant. 

At short distances, the usual asymptotically free Coulomb form (transforming as the 
zeroth component of a Lorentz 4-vector) is obtained: 

W = (20) 
3 r 

where a s (r) is obtained in the leading log approximation and is parametrized as follows 

2tt 

(11-^44 + f] 



= m^Jirm (21) 



The parameter A defines the "long distance" saturation value for a s . We use A = 
2 which corresponds to a s (r = oo) = 1.0. The parameter B is related to Ajj-g by 
B = (2.23Ajgg) _1 for Np = 3. We use Np = 3 fll3l throughout because for most distance 



regimes relevant to our calculation the cc and bb vacuum polarization contribution should 
be negligible. We generally use the present experimental average w 0.240 GeV \14 \ 



which corresponds to B — l.87GeV~ . As described in the results section, we do vary 
the value of B (i.e. of A.jfg) with respect to the above value for the specific purpose of 
studying the sensitivity of the shape of £(v.v') to this parameter. We find only a weak 
dependence. 

To describe the long distance behaviour we use a linear term in the potential that 
transforms as a Lorentz scalar (mass-like). Many theoretical and phenomenological 
arguments seem to favour this form [15], [16], [L7| , 

V L = !°ar , (22) 

where 7 is the usual Dirac matrix, rather than an admixture with a linearly rising th 
component of a vector term. We do our calculations with three different values for the 
string tension parameter a. The choices that we made were arrived at as follows. The 
experimental information available for the D and D s systems seems to indicate that the 
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(spin averaged) splitting between P states and S states is in both cases approximately 
0.45 GeV. We found that in order to obtain a splitting of this magnitude with our 
V c (r) + c + 7°ar as defined above, a had to be about a ~ 0.25GeV 2 . A previous study 
WM agrees with this calculation. On the other hand, this value for a is significantly 



larger than those obtained from the Regge slope data (qq systems) [[16], \\Jj. Of course, 
one cannot exclude the possibility that the 2P-1S splitting of about 0.45 GeV holds for 
the D, D s systems but is significantly smaller in the hypothetical limit Mq —>■ oo, which 
are the systems we are trying to describe in the present work. The difference would be 
due to heavy quark recoil effects in the charmed systems. In fact, preliminary evidence 
was reported recently for a B**(£ = 1) candidate with mass Mb** ~ 5610Mel / (Ref. |18| 



and also talk by V. Liith at the same conference), which would indicate a 2P-1S splitting 
of only about 0.34 GeV. If confirmed, this would allow us to use more conventional values 
for cr, as extracted from the Regge slope a'. We therefore also consider this possibility. 
The extraction of a from a' is somewhat model dependent. Two common relations are 

a = — !— (string model [19]) (23) 
2ixa' 

and 

cr = — (2 — body generalization of Klein Gordon Eq. [16, |T^]) (24) 
8a 



Using a' ~ 0.9GeV 2 [T!| we obtain respectively cr « 0.18GeV 2 (string model) and 
a « O.lAGeV 2 (2-body K.G. Eq.). The latter value agrees also with a recent lattice 



estimate [20|. We would like to remark that if the 2P-1S splitting for the B mesons 



(and therefore also in the mq —> oo limit) turns out to be similar in magnitude to the 



observed splitting in the D, D s systems, as predicted in several models plj we could still 
fit this splitting using the more conventional values of cr (0.18 and 0.14 GeV 2 ), provided 
that we change the parameter in V c (r) (See equation (21) and discussion below) to 
Ag)_ > O(0.5GeV) , i.e. B < O(0. VGeV- 1 ). 

So far we have used as input the well known leading short distance and long distance 
behaviour of the potential. There is less theoretical knowledge about the shape of the 
potential in the intermediate region. We introduce an additive constant term c D , which 
is clearly subleading both in the short and long distance regimes. Because the only role 
of this constant is to define the absolute scale of the light quark energy e q , this constant 
gets absorbed once we identify e q = A q = lim mQ ^oo(M g Q — itiq) and assign A q some 
physical value. 

We obtain a plausible range for the physical parameter A q from previous theoretical 
works that estimated the value of the "pole mass" of the b quark, m&. We use the 
relation 

< P 2 > 

M B (Spin averaged) w 5310MeV w A u d + m b H q - — , (25) 

2m b 
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where < Pq >~< P g 2 > has been used. 

A recent lattice study finds m& = 4950±150MeV H, which also agrees with a HQET 
estimate fl§. QCD sum rule estimates are lower, closer to m;, « 4.6Gel / f|]. From these 
values for m;, and eq (25) we then find the range 

A M » 150 - GOOMeF . (26) 

We do our calculations using mostly the two extreme values as well as the average 
value in this range = 150, 375 and QOOMeV. Notice that we have to make a 
distinction between K u>d and A s , but this does not introduce any complications because 
experimentally 

M Ds — Mo ~ M Ba -M B ^ lOOMeV , (27) 
Mn. - M Ds « Md* - M D , (28) 



and 

Thus, simply taking 



M B . - M Bs « Mb* - M B . (29) 



A s = K 4 + lOOMeV (30) 



seems to be consistent with the heavy quark expansion to 0(—). We will adopt this 
prescription to compare results for mesons containing a "u" or "d" quark with those 
containing an s quark. 

In summary, our Dirac equation potential is of the form 

V = ~ M * 7l r v=* + c ° + ^ 0(jr ( 31 ) 

27rln(2.0+ L87(G r eF } ) 

where we assumed af 3 = 1, Np = 3, Ajjg = 0.240 ; a takes the values 0.25, 0.18 or 
0.14 and c D takes values such that e u ^ = A Ujd takes a value we select in the likely range 
150-600 MeV. For the quark masses we use m u = = 0, m s = 0.175GeV. We have 
checked that this choice for m s is consistent with eq (30). 



4 Quantitative results and discussion 

In this section we use the parametrization for the Dirac equation discussed in the previous 
section in order to find the light quark wavefunction, ip q (x). We then use the formalism 
developed in section II, to calculate the Isgur-Wise function £(v.v') in the Breit reference 
frame, for different values of the parameters in the potential and for systems where the 
light quark is u,d [ra u = m d w 0) or s (m s ~ 0.175Ge\^). 
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In the case of a central potential, the time-independent Dirac equation for a state 
with angular momentum quantum numbers j and m, is reduced to radial equations by 
writing 

^ = -f 9ir l?$ 6 '*l (32) 

where k — — £ — 1 for j = £ + 1/2 and k — £ for j = £ — 1/2 and f2 Km is a spinor with 
spin 1/2 coupled to angular momentum. The radial equations are then 

[VJ(r) + m]«?(r) + [~ + £]/(r) = eg(r) , (33) 

+ -]g(r) + [-m + V.(r)]/(r) = ef(r) . 
ar r 

In the present calculation the potentials operating on the large and small components 
are parametrized as 

Vi(r) = ar - 27rln[ j r +.B/r] + °o > ,„.s 

W = -^- 27rln[W] +C °- 

The eigenvalue problem has been solved in two independent ways. The algebraic 
approach is to expand g(r) and f(r) in "basis" functions 0«(r) , i = 1,...,M and 
Xj( r ) j 3 ' = 1) leading to a matrix eigenvalue problem of dimension (M+N) for e. 
The smallest N eigenvalues correspond to hole states and the (N + l) st eigenvalue is the 
lowest particle state. In the present calculation <p and x are taken to be 



0.(r) = ^e-^ 2 / 2 , i = 1,...,M , 
Xj (r) =r j e -°r 2 /2 ) J = 1,..., AT . 



(35) 



This form is chosen to have the correct analytic behaviour at large r values. Matrix 
elements of the kinetic energy operators and the linear potential can be evaluated ana- 
lytically, but matrix elements of the "Coulomb" term are computed numerically. Satis- 
factory results are obtained for modest values of M and N of order 10. 

The problem has also been solved fully numerically in the finite difference approxi- 
mation by converting the two first order equations to the second order Pauli equation. 
Details of the method used can be found in Ref. [22]. Energies found by the two meth- 
ods are identical. However, for evaluating the form factor the more accurate numerical 
wavefunctions have been employed. 

We present our results for £ u ,d,s( v - v ') i n graphical form for the range 1 < v.v' < 
4 which is of phenomenological interest for decays of the types B, B* — > D,D* + 
X ■ B,B* -> K* + X ; D,D* -> K* + X. Figures (1), (2) and (3) exhibit our re- 
sults for "inertia" parameter A u4 = 0.6 GeV, 0.375 GeV and 0T5 GeV (A s = 0.7 GeV, 
0.475 GeV and 0.25 GeV ) respectively. For each value of A Ut d we show the result 
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iuAiv-v') for a = 0.25GeV 2 , 0.18GeV 2 and O.UGeV 2 while for each A s we give £ s (v.v') 
with a = 0.18GeV 2 for comparison with £, u ,d(v.v'). 

In Table 1 we present the values of the zero recoil slope , £,'(v.v')\ v . v ' = i for the same 
values of the parameters A U)d , A s and a as used in figures (1-3). 

We observe that for A u ^ ps 0.15, £, u ,d( v - v ') is almost independent of the parameter a. 
It is clear from eqs (16-18) that if |A|(= |e 9 |) is very small, £(v.v') is controlled mostly by 
purely kinematic factors, the shape of the wavefunction becomes unimportant (provided 
that it is properly normalized). We observe also that £ u ,d(v-v') and £ s (i>.i>') are quite 
close (a = 0.18 , A s PS A U)d + lOOMeV) for all values of A Ujd . This confirms that the 
strange quark can be treated as a light quark and that SU (3) p is only softly broken for 
the processes considered here, once the (phenomenologically imposed) shift in the value 
of A has been considered. 

Although it is not obvious from figures (1-3), we checked (by using different values 
of a and A U;d and keeping A 2 ld /a constant) that £ u ,d( v - v> ) m the range 1 < v.v' < 4 
depends on the ratio A 2 ld /a only, to a precision better than 2%. This "scaling" is not 
completely unexpected since a and A Utd are the main dimensionful parameters entering 
the problem, and £(v.v') is dimensionless. The other dimensionful parameter is Ajfg 
in V c , which we are keeping constant, affects V c only logarithmically and, for massless 
quarks, does not much affect the shape of the wavefunction far from the origin (Recall 
that a pure vector-like potential does not produce bound states for massless quarks, and 
we are using m u ps m d ps 0). 

Figure 4 illustrates this scaling effect by giving the values of the zero recoil slope 

A 2 

p 2 = — £'(1) for different values of The points were obtained from Table 1 using 
the 3 values of a and A u>d quoted. It turns out that all the points satisfy to a very good 
approximation the empirical linear relation 

P 2 = -C,,(l)^ + 0.39^ (36) 

A 2 

where the numerical coefficient of -^varies only by ±0.01 within our range of A Ujd and 
a. This equation can then be used to find the zero-recoil slope within our model for an 
arbitrary value of the imput parameters A u>d and a. 

We checked whether this approximate dependence of £(v.v') on the ratio A 2 d /<j only 
also holds for the system with a strange (light) quark (m s ps 0.175GeV). We found 
that C, s (v.v') does depend on a and A s independently, i.e. not only through the ratio 
A 2 /a. What causes the different behaviour in this case is that there is an additional 
dimensionful parameter in the problem, m s , so that we can form two independent di- 
mensionless ratios A/m s and c/m 2 - Regarding p 2 = —^(1), we notice from Table 1 that 
consistently p 2 > p 2 u d , i.e. the Isgur-Wise function has (exhibits) a more rapid decrease 
for qQ mesons where the light quark is an s quark. This result was obtained as well in 
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the calculations of Refs f7|, |23 . 



We also studied the effects of changing the value of Ajj$ (B parameter) in V c (see 
eq 20, 21 and discussion below) on £(i>.i/), keeping a and A fixed. We found that 
the behaviour of £(v.v') is not very sensitive to changes in Aj^. For instance, using 
a = 0.18GeV 2 and A M = 0.375GW (A, = 0.475GeK), a change of A^ by 50% either 
way from our central value A^g^ « 0.2AGeV changes the zero recoil slope £'(1) by only 
about 2.2% (3%). The relative changes in (1— £(i>.t/)) are of the same order of magnitude 
for the whole range 1 < v.v' < 4. 

Finally, we compare our results for £ u ,d( v - v ') with the existing experimental data and 
other theoretical estimates. 

As mentioned above, we observed that ^u,d( v - v ') m our model depends (to a very 
good approximation) only on the ratio A 2 ld /a. Hence, we adopt a fitting procedure to 



the ARGUS -93§ and CLEO -93|Tg and CLEO '94 |TJ data in which the parameters 
of the potential, A=2, B=1.87 GeV~ x and a = 0.18GeV 2 are kept fixed, while A u ^ and 
the physical observable |V^,|/i are varied in such a way that \ 2 i defined below as being 
proportional to the total square deviation from the data (weighed by the experimental 
uncertainty for each point) is minimized: 



N 



X 2 = E (\VcMKv.vV - fiftf , (37) 

where 



H = J Y^- s ( nt s to ARGUS '93 and CLEO '93) (38) 

and 

fi = F(l)^riA in the notation of [11] (fits to CLEO '94) (39) 

The normalization factor 1/(JV — 2) (N experimental points, 2 fitting parameters) allows 
us to interpret the magnitude of x 2 f° r the best fit, Xo, as an absolute measure of the 
"quality"of the fit to the corresponding data set. The factor F(l) ~ tja required to 



extract an estimate of | V c t\ from the CLEO '94 data [IT] has been estimated theoretically 
by several authors: 



F(l) = 0.97 ±0.04 [1, 
F(l) = 0.96 ±0.03 |1 (40) 



F(l) < 0.94 , F(l) w 0.89 ±0.03 p6 



For each data set, we obtain "acceptable" ranges for our two fitting parameters, A n ^ 
and \V c b\n by requiring 

X 2 < J^xl , (41) 
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where Xo corresponds to the best fit (i.e. Xo is the minimum of the function 
X 2 (Au,d , Because of the scaling behaviour described above eq (36), the ranges 

that we obtain for the parameter A U)d and fixed a = 0.18 GeV 2 can be translated into 
ranges for the dimensionless ratio J a where a is allowed to take values other than 
0.18 GeV 2 . 

We would like to add a cautionary note before giving the results of our fits. Because 
the Isgur-Wise function that we calculate in our model and use for these fits does not 
take into account finite-MQ corrections (0(Aq CD /Mq) at v.v' = 1 and 0(Aqcd/Mq) 
elsewhere), we have to be aware that direct comparison of our model with experiment 
(via our x 2 function) can be reliable to leading order in Aqcd/Mq only. At v.v' = 1, 
this uncertainty is of order Aq CD /Mq and can be absorbed into the parameter /i which 
multiplies \V c b\ (see eqs 38-40). However, away from v.v' = 1 there is an "intrinsic" 
uncertainty in the comparison of our model to experiment of expected relative magnitude 
of order (v.v' - \)A QC d/Mq. 

The best fit to the ARGUS '93 data (8 points) Q gives xl = 0.54. Using the 
above prescription (eq (41)) we obtain the following acceptable ranges for the fitting 

A 2 

parameters: = 4.8 ±1.7 (corresponding to (see eq (36)) p 2 d = —£,' u d (l) = 2.4 ± 0.7) 
and = 0.050 ± 0.008. We note that even if a small value for the string 

tension a = O.lAGeV 2 were used, the acceptable range for the inertia parameter would 
be A Ui d ~ 0.81 ±0.15 GeV, which is significantly above most theoretical estimates (see eq 
(26)). The corresponding range for p 2 ud overlaps with some of the theoretical estimates 
(see Table 2) but is centered above most of the predicted ranges. 

The best fit to the CLEO '93 data (7 points) Jl0[ is poorer with xl ~ 1-18. Here we 

A 2 

find the ranges = 3.4±1.8 (corresponding to (see eq (36)) p 2 ud = — £' u d (l) = 1.8±0.7) 
and iKbly^If^ = 0-043 ± 0.008. There is some overlap between the resulting range for 
A u>d (eg A Ujd = 0.66 ± 0.19 for a = 0.14 GeV 2 ) and independent theoretical estimates 
(See eq 26). Also, the range for the slope p 2 u d significantly overlaps with several previous 
theoretical predictions (See Table 2). We would like to remark that if we ignore the 
CLEO '93 data point corresponding to highest recoil (v.v' « 1.5), the fit to the remaining 

A 2 

6 points is greatly improved. We obtain in this case xl — 0.62 and ranges = 2.1 ±1.6 
(corresponding to p 2 ud = 1.3 ±0.6) and iKftlyj^j = 0.038 ±0.006. For commonly used 

values for the string tension (a = 0.14 — 0.18 GeV 2 ) the acceptable range for A Ujd in 
this case is centered well within the range of previous theoretical estimates (See eq 26). 
The range for the zero recoil slope p 2 d = —£' u d (l) is very similar to the ranges obtained 
in recent lattice estimates |3T|, [32|] as well as from several other theoretical calculations 
(See Table 2). 

The best fit of our model to the recent CLEO '94 data analysis (7 points) (Tl| gives 
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Xq = 0.50 which is the lowest Xo °f ai l our fits, in spite of the smaller experimental error 
bars. This means that within a region of its parameter space, our model agrees well 
with this data set. In turn, we expect that this data with smaller error bars will be 
useful in selecting a relatively narrow region for our model parameters (mainly J a) 
as well as for the Standard Model flavour mixing parameter \ V c b\ (or at least the product 

A 2 

|V^|-F(1)). The parameter ranges obtained by imposing eq (41) are: — ^ = 1.9 ± 0.7 
(corresponding to p 2 ud = -& >d (l) = 1.3 ± 0.3) and |K fe |F(l) = 0.037 ± 0.003. As 
expected, the smaller experimental error bars lead to a better determination of our 
model parameter A 2 d /a (and correspondingly, via eq (36), of the slope at zero recoil 
p 2 d ) as well as of the physically interesting product |V^b|F(l). We note that for the 
commonly used values for the string tension a = 0.14 — 0.18 GeV 2 , the resulting range 
for A U) d (0.41 GeV < A U)d < 0.69 GeV) overlaps significantly with the upper half of 
the range given in eq (26), which was obtained from recent theoretical estimates of mj 
0, H |]. In this respect, the larger values of a (eg o = 0.25 GeV 2 (See discussion below 
eq (22)), seem to be less favoured. We would like to remark here that an estimate of 
the pseudoscalar decay constants fr> and Jb in the context of the relativistic model used 
here [j39H , when compared with recent estimates of these constants with lattice and QCD 
Sum Rule methods, favours the lower values for a as well. 

Because the Isgur-Wise function is normalized at zero recoil, £(1) = 1, the slope at 
zero recoil £'(1) determines to a good approximation the value of the function close to 
v.v' = 1. Therefore, at least close to v.v' = 1 (all the existing data is in the interval 
[1,1.5]) the slope £'(1) = — p 2 is a reliable tool for comparison of the different theoretical 
estimates of £(v.v'). Our original model parameter ranges 0.15 GeV < A u ^ d < 0.6 GeV 
and 0.14 GeV 2 < o < 0.25 GeV 2 (See discussion between eqs (22) and (26)) give 
(through our result in eq (36)) a wide range 0.54 < p 2 d < 1.5, which overlaps with 
many different theoretical estimates (See Table 2 and refs [6-8] and [27-38]). The only 
exceptions are refs [33-36]. On the other hand, the relatively narrow range for p 2 ud 
favoured by our fit to the CLEO '94 data (p 2 ud = 1.3 ± 0.3), overlaps with the 
predictions in refs [7], [8], and [28-32] only, although is not far from the range advocated 
in ref [6] (See Table 2). Both the lattice gauge theory calculations [31] and [32] contain 
the range favoured by the fit of our model to the CLEO '94 data within their predicted 
ranges. 

The extraction of a precise value for the Standard Model parameter \V c b\ from the 
range \V cb \F(l) = 0.037 ± 0.003 favoured by our fit to the CLEO '94 data is partially 
hindered by the present uncertainty in the theoretical determination of F(l). For ex- 
ample, if we use the whole range of values for F(l) given in eq (40), we would obtain 
0.034 < |V c6 | < 0.046. 

It is interesting to note that our fits to the CLEO '94 data (7 points) produce similar 
central values for A 2 ^ d /cr, p 2 ud and |V C {,| (for F(l) < 1) as do our fits to the CLEO '93 
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data with the point corresponding to the largest v.v' omitted from the CLEO '93 set (i.e. 
6 points). The resulting ranges for these quantities are, however, significantly narrower 
(by a factor of about 2) for the CLEO '94 data set. We should also point out that the 
above mentioned central values for (? ud and |V c f,|F(l) that we obtained from the fits of 
our model to the CLEO '94 data set (Fig 12 in ref fill ) are both somewhat larger than 
the central values for these quantities obtained in the CLEO analysis carried out in ref 
UTTf (In their notation a 2 = p^ d ). Once the uncertainties are included, however, our 
results for these quantities are compatible. 

Eventually, when the experimental data becomes more precise and the 1/Mq correc- 
tions can be incorporated with fewer uncertainties, one should be able to further narrow 
the allowed range for our main model input, A u ^/a, as well as for the flavour mixing 
parameter \V c b\. 
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Table Captions 

Table 1 : The slope of the Isgur-Wise function at zero recoil for various values of 
the parameters A and a. 

Table 2: Various theoretical estimates of p 2 u d = —£'(1). 



20 



Figure Captions 

Figure 1: The Isgur-Wise function f M (&) for A M = 0.15 GeV (A s = 0.25 GeV) 
and cr = 0.25 , 0.18 and 0.14 GeV 2 (a = 0.18 GeV 2 ). 

Figure 2: The Isgur-Wise function f U)d (&) for A M = 0.375 GeV (A s = 0.475 GeV) 
and a = 0.25 , 0.18 and 0.14 GeV 2 (a = 0.18 GeV 2 ). 



Figure 3: The Isgur-Wise function f M (&) for A M = 0.6 GeV (A s = 0.7 GeV) 
and o- = 0.25 , 0.18 and 0.14 GeV 2 (a = 0.18 GeV 2 ). 

Figure 4: The slope of the Isgur-Wise function at zero recoil p 2 u d = — £„ d (l) for 

A 2 

different values of the parameter 



21 



Table ' 






a = 0.25 


a = 0.18 


a = 0.14 


Am = 0.15 

A s = 0.25 


P'U = 0-53 
p 2 s = 0.58 


= °- 55 

= 0.61 


Pl, d = 0-56 
= 0.64 


K,d = 0-375 
A s = 0.475 


P'U = 0-72 
pi = 0.80 


PU = 0-81 
= 0.90 


Pl, d = 0-89 
= 0.99 


K 4 = 0.6 
A s = 0.7 


Pl A = 1-06 
P 2 S = 1-16 


PU = 1-28 
p2 = 1.37 


pIa = i- 50 

pi = 1.56 



Table 2 



Bjorken [8] 


>i 


Isgur et al. [27] 


0.63(0.33) 


Rosner [28] 


1.44 ±0.41 


Mannel et al. [29] 


1.77 ±0.74 


Neubert [30] 


1.28 ±0.25 


Bernard et al. [31] 


1.41 ±0.19 ±0.41 


UKQCD Collaboration [32] 


1.23 


Radyushkin [33] 


oo 


Karanikas and Ktorides [34] 





Sadzikowski and Zalewski [7] 


1.24 


Kugo et al. [35] 


1.8-2.0 


Ivanov et al. [36] 


0.43 


Ivanov and Mizutani [37] 


0.42 - 0.82 


Narison [6] 


0.52-0.92 


Blok and Shifman [38] 


0.5-0.8 
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